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Abstract: 


The process of heat exchange in mechanical systems is discussed in the article. The solution to two- 
dimensional non-stationary heat conduction problems in the presence of point heat sources is 
studied based on the finite element method. To solve a plane heat conduction problem, a triangular 
finite element with rectilinear sides is used. The basic relations for linear triangular finite elements 
are presented, taking into account the influence of point heat sources. A source set at a specific 
point on the plane describes a heater in the form of a string. This point on the plane is the trace of 
the string and is set by the heat generation power per unit length. Based on the developed 
computational algorithm and software, a numerical solution to a specific problem is given and the 
influence of point heat sources on the distribution of the temperature field of the structure is 
studied. 
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In many engineering problems, knowledge of the temperature distribution in a structure is an 
important aspect. The amount of heat supplied to it or lost by it can be calculated provided that the 
temperature distribution is known. In the article, the temperature field is formed by considering the 
influence of a point heat source set at a specific point on the plane and describes the presence of a 
heater in the form of a string. The trace of the string is a point on the plane, determined by the 
power of heat generation. 


Formulation of the problem 


The problem of temperature distribution at different time points is solved in a plane formulation [1, 
2]. To find the temperature field in the rectangular calculation domain (a plate) ABCD (Fig. 1), a 
two-dimensional non-stationary heat conduction problem is solved based on equation [4]: 


2, 2. 
ey ae ly qe ie) (1) 
Ot Ox ~ Oy 
where A=cp is the specific volumetric heat capacity; c is the specific heat capacity of the 
material; p is the density; K,.,K,, are the heat conductivity coefficients in the corresponding 


directions; Q is the power of the heat source inside the body. 
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Fig. 1. Scheme of the computational domain 
Symmetry conditions are set on lines AD and BC (the so-called natural boundary conditions [3]): 
OT or 
=(), 


al, a 


=); (2) 


BC 


A constant temperature (Dirichlet condition) is set on line AB, equal to the ambient temperature T,,: 


Aes (3) 
Conditions that determine the convective exchange of heat with the medium can be set on line AB: 
Ky, — she 1) =0, (4) 

oy 
where / is the heat transfer coefficient. 
The temperature field throughout the entire computational domain at time t = 0: 
T(x, y,0) =T°(x, y) is set as the initial condition. 


In this case, the solution to the system consisting of equation (1) and boundary conditions (2) - (4) 
is reduced to minimizing the functional [4]: 


1 ery or) eT h 
=|2 Kh. |=) 2k.) —— | =807 324-2? lav +1 0 = 5 
x FI (=) (5) QO a | 5° vy, © 


where S, is the area of surface where convective heat exchange occurs. 
Finite element solution to the problem 
The temperature in the finite element is specified as the product of two independent functions: 
T° =N(x, y)TO) 
or in a matrix form: 
T(t) 
T =| N,(x,y)N On y)N, 0% y) 47,0 p (6) 
T(t) 


For finite element e, the condition for the extremum of the functional leads to the following 
system of differential equations: 
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fe eer} fr =o. a) 

here 

[el= | ANILNT av, (8) 
‘. 

[k°] = [ [BLD IBY av + | ALNILN Tas, (9) 
ve So 

(F}=— | QAN°TdV— J he, [N° Tds, (10) 

yee) so 
where 


v¢ is the volume of the finite element; 

[N] is the shape function matrix; 

[B*] is the matrix of shape function derivatives; 

[D*]is the matrix of material properties containing thermal conductivity coefficients. 

Integral (8) when solving a two-dimensional heat conduction problem has the following form 


(damping matrix): 


_ AAa, 


. 1 
[c*] ) 
1 


I 
2 1), (11) 
iL 2 


where 
A is the area of the finite element; a, is the element thickness. 


The thermal conductivity matrix (in the absence of convection) has the following form: 


‘ bb, bb, bb, Cl, CC; CC, 
a 
[ke] =n bb, bib, bb, \+ re CO GE, OG. |: (12) 
bb, bb, bb, Cyl, YC; Cy Cy 


where b =y,—y,3 ¢; =X; —%,5 (the remaining quantities are obtained by circular permutation of 
indices 1, j, k). 


If side i of the j finite element is subject to convection, then the second integral in (9) is: 


210 
T ha,L, 

[ALN IN] dS=—“|1 2 oO], (13) 

82 00 0 


where L, is the length of the side of the element between nodes i and j. 
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If heat is lost by convection between sides with nodes j, k or k, i, then matrix (13) is transformed 
[4]. 


We assume that the value of Q is constant inside the element, then: 


ie 1 
QfINIav =QffL, Jav= S21 (14) 
Vv Vv L , 1 


3 


If convective heat exchange occurs on the side with nodes i j, then the “load vector” in solving the 
heat conduction problem is: 


a ATL, 
[AT.1N] dS =—*4 


0 
For a finite element mesh, a system of ordinary differential equations is written as: 
cS KUT HLF} =0, (15) 
where 


[l= > lel (kis > 1k); T=, 


Replacing the time derivative in equation (15) with its finite-difference analog, we obtain an 
implicit difference scheme for solving the heat equation using the finite element method [3-6]: 


[fel txa)iry = (e (KITV ~(Fy"" (16) 

Thus, if the temperature vector {r\" at time point ¢, is known, then the temperature of the plate at 
time point ¢,,,=f,+Ar is obtained as a result of solving the system of linear algebraic equations 
(16). 

Computational experiment and analysis of results 

The following problem is considered as a test example [7]: 


Three cables pass through a heat-conducting medium, as shown in Fig. 2. The thermal conductivity 
W 


cm: K 


coefficients of the medium is K, =K y= 50 . The heat transfer coefficient on the surface of 


the medium is h=10 . The medium under consideration is limited on the sides by a thick 


cm’ - K 
layer of insulation. Air temperature on the surface of the medium is 7, =+40°C. Temperature of 
the lower layer of the medium is 7, = +20°C. The heat radiation power of each cable is Q = 200W 
Copper has the following § thermophysical characteristics (1, in Fig. 2): 
A= 384 W/(m°C), p = 8800 kg/m’ ,c= 381 J/(kg°C) . Thermophysical parameters of the material of 
inclusion are: steel (2, in Fig. 2): A, =46W/(m-°C), p, =7800kg / m’°, cy = 460J / (kg -°C). 
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Dimensions of the material of inclusion: n=18 cm, m=2.375 cm. Plate dimensions: a=5 cm, b=4 
cm, c=5 cm, d=5 cm, e=4 cm. 


Fig. 2. Arrangement of cables in a heat-conducting medium 


The temperature distribution in a homogeneous material (Fig. 2.a) at t = 6, 120, and 300 seconds is 
presented in Fig. 3. At the initial stage of the process (t = 6 sec), there is no influence of sources, 
then at (t = 120 and 300 sec), the influence of sources is observed, characterized by thermal 
concentration in the vicinity of sources [9-12]. 


C) 10,12 20 

1,67 11,21 20,57 
3,33 12,31 21,14 
5 13,4 21,71 
6,67 145 22,28 
8,33 15,59 22,85 
10 16,69 23,42 
11,67 17,79 23,99 
13,33 18,88 24,55 
15 19,98 25,12 
16,67 21,07 25,69 
18,33 22,17 26,26 


t=6 sec t = 120 sec t = 300 sec 


Fig.3. Temperature distribution in a homogeneous material 


In the presence of a rectangular steel area in a structure made of copper (Fig. 2 b), a redistribution 
of the generated heat is observed in its lower part (Fig. 4). The change in the temperature field 
reflects the difference in the thermal conductivity coefficients of copper and steel. 


t =6 sec t =120 sec t =300 sec 


Fig.4. Temperature distribution in a non-homogeneous material 
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For a specific point (2.5 cm, 11.5 cm) located on the left side of the cavity, table 1 presents 
temperature values that indicate the influence of the non-homogeneity of the structure, namely a 
decrease in temperature in the vicinity of sources. 


Table 1. Temperature at point (2.5 cm, 11.5 cm) 


Plate material t=6 sec t=120 sec t=300 sec 

Homogeneous 0.0358 °C 17.6320 °C 30.3977 °C 

Non-homogeneous 0.0016 °C 10.5583 °C Za02)0 Ek 
Conclusion 


Thus, the presence of sources increases the temperature field in the vicinity of sources, and the 
presence of non-homogeneity changes the overall temperature background of the structure. 
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